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ABSTRACT 

The z = 2.286 IRAS galaxy F10214+4724 remains one of the most lumi- 
nous galaxies in the Universe, despite its gravitational lens magnification. 
We present optical and near-infrared spectra of F10214+4724, with clear 
evidence for three distinct components: lines of width ~ 1000 kms -1 from 
a Seyfert-II nucleus; 5; 200 km s _1 lines which are likely to be associated 
with star formation; and a broad (~ 4000 km s^ 1 ) Cm] 1909A emission 
line which is blue-shifted by ~ 1000 km s^ 1 with respect to the Seyfert-II 
lines. Our study of the Seyfert-II component leads to several new results, 
including: (i) From the double-peaked structure in the Lya line, and the 
lack of Ly/3, we argue that the Lya photons have emerged through a neu- 
tral column of Ah ~ 2.5 x 10 25 m~ 2 , possibly located within the AGN 
narrow-line region as argued in several high redshift radiogalaxies. (ii) The 
resonant Ovi 1032, 1036A doublet (previously identified as Ly/3) is in an 
optically thick (1:1) ratio. At face value this implies an an extreme density 
(n e ~ 10 17 m~ 3 ) more typical of broad line region clouds. However, we 
attribute this instead to the damping wings of Ly/3 from the resonant ab- 
sorption, (iii) A tentative detection of Hell 1086 suggests little extinction 
in the rest-frame ultraviolet. 

Key words: galaxies: active - galaxies: formation - galaxies: individual 
(FSC 10214+4724) - galaxies: starburst - gravitational lensing 



1 INTRODUCTION 

The z = 2.286 IRAS galaxy FSC 10214+4724 is one of 
the most apparently luminous objects in the Universe, 
and its discovery (Rowan-Robinson et al. 1991) led to 
much speculation about its possible status as a pro- 
togalaxy. This speculation was based on the extreme 
bolometric luminosity of the object (Rowan- Robinson 
et al. 1991), and more specifically on the huge gas 
mass and star formation rate inferred from the sub- 
mm molecular line and continuum detections (e.g. 
Solomon, Downes & Radford 1992; Rowan- Robinson 
et al. 1993). 

This speculation was dampened by a series of pa- 



pers which proved that F10214+4724 is being grav- 
itationally lensed, and, at all wavebands, is intrinsi- 
cally an order of magnitude or more dimmer than 
it first appeared. Although a lensing bias was sus- 
pected by a number of authors (e.g. Elston et al. 1994, 
Trentham 1995), the first direct empirical evidence of 
strong gravitational lensing was provided by a deep 
near-infrared image (Matthews et al. 1994) which re- 
vealed an arc-like structure centred on a galaxy close 
to the line of sight to F10214+4724. Several sets of 
authors published lensing interpretations (Graham & 
Liu 1995; Serjeant et al. 1995; Broadhurst & Lehar 
1995) which were confirmed by the appearance of the 
HST image of Eisenhardt et al. (1996): this image con- 
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tained highly elliptical, high surface brightness fea- 
tures characteristic of strong lensing, as well as a clear 
counter-image. These papers also attempted to con- 
strain the redshift of the system responsible for the 
gravitational lensing, our contribution (Serjeant et al. 
1995) being spectroscopy of two galaxies projected 
w 1 and ^3 arcsec from F10214+4724. This work 
revealed tentative 4000A breaks at z ~ 0.90 in both 
galaxies, later confirmed using fundamental plane ar- 
guments in HST imaging (Eisenhardt et al. 1996), 
and also tentatively supported by a weak absorption 
line at z = 0.893 in the F 10214+4724 spectroscopy 
of Goodrich et al. (1996). The HST R-band image of 
F10214+4724 implies magnifications of ~ 100 (Eisen- 
hardt et al. 1996) in this waveband, but it now seems 
likely that differential flux magnification causes lower 
magnification factors for the more extended structure, 
as argued by several authors. 

Gravitational lensing appeared to offer a com- 
pelling explanation for the extreme luminosity of FSC 
10214+4724 (e.g., Broadhurst & Lehar 1995). As a 
result, the IRAS galaxy is no longer so extreme in 
its properties: indeed, in many respects, it resembles 
local ultraluminous infrared galaxies and Seyfert-II 
galaxies. Nevertheless, both Downes et al. (1995) and 
Green & Rowan- Robinson (1996) argue for a bolo- 
metric magnification factor of < 10 using arguments 
based on minimum black body sizes. FSC 10214+4724 
remains one of the most intrinsically luminous objects 
in the Universe. The importance in this object still 
lies in studying whether high-z hyperluminous activ- 
ity, such as that seen in FSC 10214+4724, differs in 
all but scale from local objects, and in determining 
the relative contributions of the starburst and AGN 
components (Elston et al. 1994; Lawrence et al. 1993, 
1994; Soifer et al. 1995; Goodrich et al. 1996; Kroker 
et al. 1996; Hughes, Dunlop & Rawlings 1997). 

Prior to making the observations reported in 
this paper there had been no direct evidence for the 
presence of an embedded broad-line (e.g. Seyfert- 
I or quasar) nucleus in F10214 + 4724, although 
high («20 per cent) rest-frame ultraviolet polariza- 
tion (Lawrence et al. 1993; Jannuzi et al. 1995) sug- 
gested that one was present. This situation changed 
with the deep spectropolarimery of Goodrich et al. 
(1996) showing clear broad lines in polarized light. 
This extended the close spectral similarities between 
F10214 + 4724 and Seyfert-II galaxies, specifically 
NGC1068, which was first remarked on by Elston et 
al. (1994). 

Prior to our observations there had also been 
no reported detection of the narrow (;$ 200 km s _1 ) 
emission lines expected from any star-forming activ- 
ity in F10214 + 4724. A star-forming component is 
expected if the analogy with NGC1068 is to be com- 
plete. This situation also changed during the prepara- 
tion of this paper. Using a novel imaging near-infrared 
spectrometer Kroker et al. (1996) presented evidence 
for spatially-extended narrow Ha emission just as ex- 



pected if the Seyfert-II nucleus of F10214+4724 is ac- 
companied by a circumnuclear starburst. 

In this paper we present, analyse and interpret 
optical and near-infrared spectroscopy of IRAS FSC 
10214+4724. The details of data acquisition and anal- 
ysis are given in Section 2. In Section 3 we compare 
our results with previous and contemperaneous spec- 
troscopic studies of F10214+4724. In Section 4 we 
interpret the data on the Seyfert-II emission line re- 
gion including a discussion of optical depth effects on 
the resonance lines, and some modelling of the spec- 
trum using the photoionisation code CLOUDY (e.g. 
Ferland 1993, 1996). In this section we reach conclu- 
sions about the Seyfert-II properties of F10214+4724 
which differ significantly from those reached by previ- 
ous studies. In Section 5 we interpret the data on the 
region responsible for the narrow (< 200 kms^ 1 ) Ha 
line seen in our near-infrared spectrum: this feature is 
likely to be a signature of star formation. In Section 
6 we pass some concluding remarks on the nature of 
F10214+4724. Further data on galaxies foreground to 
F10214 + 4724 are given in Appendix A, and our at- 
tempts at identifying the weak emission line at 2067A 
are discussed in Appendix B. 

2 DATA ACQUISITION AND ANALYSIS 

2.1 Optical spectroscopy 

We observed F10214+4724 on the nights 1995 Jan- 
uary 28 and 1995 January 30 with the ISIS spectro- 
graph on the William Herschel Telescope (WHT) tak- 
ing advantage of 0.7-arcsec seeing. We refer the reader 
to Serjeant et al. (1995) for full details of these ob- 
servations. The spectrophotometry is in good agree- 
ment with the continuum measurements tabulated by 
Rowan-Robinson et al. (1993) as well as (accounting 
for slit losses) the line parameters of Goodrich et al. 
(1995). There are significant disagreements between 
our measurements and those tabulated by Elston et 
al. (1994) — their line fluxes are typically larger by a 
factor ~ 2 — but these can plausibly be attributed to 
imperfect correction for the non-photometric observ- 
ing conditions experienced by Elston et al. Neither 
the continuum, nor any of the emission lines are re- 
solved spatially (with « 0.35 arcsec pixels); limits on 
the extended line emission are discussed by Serjeant 
et al. (1995). 

2.2 Near-infrared spectroscopy 

We observed F10214+4724 on two occasions with the 
short (150mm) camera and 3" (single pixel) slit of the 
62 x 58 InSb CGS4 array (Mountain et al. 1990) on 
the UKIRT. The slit was aligned at a position angle of 
90°. The first of these observations, on 1992 February 
5, used the 75 lines mm -1 grating to obtain a first- 
order spectrum centred near 1.6/im; a series of 30s 
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Figure 2: Cm] line modelled as a narrow line and 
narrow Sim doublet. The individual model compo- 
nents are plotted as dotted lines, and the sum as a 
full line. Clearly this makes an extremely poor fit 
to the data. The positions of the Aim emission line 
and the z = 1.361 Mgn absorber are also marked. 

exposures were arranged in sets of four shifted by 0, 
0.5, 1 and 1.5 in wavelength — this provided Nyquist 
sampling and ensured that a given wavelength was 
sampled by two pixels. The standard 'ABBA' nod- 
ding pattern with a nod of 24 arcsec (8 detector rows) 
was used. The total exposure time was 9 x 4 x 2 x 30s 
for each point in the 126-pixel spectrum. The sec- 
ond observation, on March 23 1992, employed the 150 
lines mm grating to obtain a second-order spectrum 
centred near 2.16/im. The nod in this case was 30" 
(10 rows), and the total exposure time 10 x 4 x 2 x 
40s per point. Spectra were wavelength calibrated us- 
ing Xenon/ Argon lamps and night- sky hydroxy 1 lines, 
flux calibrated using HD105601, and the positive and 
negative channels were extracted from their respec- 
tive 3-arcsec wide rows. The wavelength calibration 
is accurate to ~ 0.003/xm in H and 0.001/xm in K, 
and the spectral resolving powers were 250 at H and 
1400 at K. 



3 RESULTS AND COMPARISON WITH 
OTHER SPECTROSCOPY 

3.1 Optical spectroscopy 

Figure 1 shows our WHT optical spectrum of 
F10214+4724; line fluxes, widths and redshifts are 
tabulated in Table 1. Deeper spectra covering the re- 
gion above 3900 A have been presented by Soifer et 
al. (1995) and Goodrich et al. (1996), and an MMT 
spectrum, comparable in sensitivity and wavelength 
range to our data, is presented by Close et al. (1996); 
see also Rowan- Robinson et al. (1991) for the discov- 
ery spectrum. The velocity profiles of all the emis- 
sion lines except Lya are, with FWHMs around 1000 
km s _1 , extremely similar. These line withs are small 
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Figure 3: Cm] line modelled as two gaussians. 



compared to quasar broad lines, but large compared 
to typical AGN narrow lines which are typically a few 
hundred km s" 1 (e.g. Nelson & Whittle 1996). 

The accuracy of our wavelength calibration 
means that the identification of the bright doublet 
in the far blue with Ovi is secure, and we are forced 
to conclude that Close et al. (1996) were mistaken in 
preferring Ly/3 as the identification for this feature. 
This is probably the most clearly resolved example of 
the Ovi doublet (see e.g. Kriss et al. 1992a,b; Laor et 
al. 1994), a factor we will exploit in Section 4. Above 
3900 A comparison with the published spectra shows 
that all the labelled features are real. The identifi- 
cation of the (rest-frame) 2470 A feature, which is 
also seen in NGC 1068, is new but probably uncon- 
troversial. The identification of the line at « 1805A 
is more uncertain: we have followed Lacy & Rawlings 
(1994) by identifying it with Mgvi, a li ne p resent in 
the CLOUDY models discussed in section |4.2| . A strong 
line at this wavelength is also seen in NGC 1068 (Sni- 
jders, Netzer & Boksenberg 1986). We prefer the Mgvi 
identification to either the Sill 1814. OA multiplet or 
the [Nelll] 1814. 6A line: these were suggested as pos- 
sible identifications by Snidjers et al., and adopted 
for F10214+4724 by Soifer et al. 1995 and Close et 
al. 1996. Note that the ionisation potential of Mg 5+ 
is 141.2 eV, and thus even more extreme than the 
113.9 eV required to produce 5+ . 

This leaves one unidentified line at a rest-frame 
wavelength of 2067 A. This line is also seen in the HST 
spectrum of NGC1068 (Antonucci, Hurt & Miller 
1994) where it lies in a region confused by underlying 
Fell multiplets (see also Wills, Netzer & Wills 1980 
for a discussion of a 2080 A feature in the spectra of 
quasars), but in F10214+4724 it is clearly a distinct 
narrow feature. Our attempts at identifying this line 
are discussed in Appendix B. 

The broad asymmetric base to the Cm] 1909 line 
is shown in more detail in fig 3. This broad base inde- 
pendently present (though marginally) in both night's 
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spectra. The data are reproduced by the sum of two 
Gaussian components with a broad (~ 4000 km s _1 ) 
component blueshifted by about 1000 km s _1 rela- 
tive to a component with a similar FWHM (e.g. 
ft; 1000 kms -1 ) to the other optical lines. Although 
lines of Silll] are expected to be present at some level 
in the blue wing of Cm] 1909 (see fig 2) we were un- 
able to obtain as good a fit using a superposition of 
« 1000 kms -1 lines (fig 2). This implies that the 
previous detection of hint of a broad Cm] 1909 line 
in integrated light (Goodrich et al. 1996) is not at- 
tributable solely to the neighbouring narrow emis- 
sion lines. There is also an apparent blue wing to the 
Cll]2326 line, but this feature is sensitive to the cor- 
rection for the atmospheric A band. The lack of broad 
components to other lines may be explainable by a 
combination of scattering albedo and reddening to- 
wards the scattering surface (or between the BLR and 
the scatterers, or within the scattering region itself); 
also lower signal to noise may mask broad features in 
the Lya and Nv region (figure 4). 

The Lya line has a number of unusual properties 
which have been noted by other authors. First, it is 
weak relative to the other ultraviolet emission lines, 
e.g. Nv. This appears to reflect anomalously weak 
Lya emission: ignoring Lya, both NGC 1068 {e.g. 
Snidjers et al. 1986, Kriss et al. 1992) and the IRAS- 
detected Seyfert 2 galaxy NGC 3393 (Diaz et al. 1988) 
have similar ultraviolet spectra to F10214+4724. Sec- 
ondly the Lya line is double peaked; this has been 
interpreted as requiring 'self-absorption' (Soifer et al. 
1995; Close et al. 1996). We will discuss these points 
further in Section 4. 

The continuum depression shortward of Lya is 
within the range, though at the upper limit observed 
in quasars at similar redshifts {e.g. Warren et al. 
1994), consistent with the expected lack of a QSO 
proximity effect {e.g. Giallongo et al. 1996 and refs. 
therein). The Hell 1086 line may be marginally de- 
tected. We confirm the presence of a Mg II 2800 ab- 
sorption doublet at z — 1.316 (see fig 2), but have in- 
sufficient sensitivity to confirm the weaker z = 0.892 
doublet tentatively identified by Goodrich et al. 
(1996). 



3.2 Near-infrared spectroscopy 

Our H-band spectrum of F10214+4724 is shown in 
fig 5, and our K-band spectrum in fig 6. Other near- 
infrared spectroscopy of F10214+4724 has been pub- 
lished by Elston et al. (1994), Soifer et al. (1995), 
Iwamuro et al. (1995) and Kroker et al. (1996). 

We obtain an approximate lower limit on the 
[Om]4959+5007 to H/3 ratio of 10. Iwamuro et al. 
(1995) claim a tentative detection of H(3 with an Om 
to H/3 ratio of 27, though the H/3 is a factor of 3 
brighter than the upper limit quoted by Elston et al. 
(1995). Iwamuro et al. used the OH airglow supres- 
sor spectrograph on the UH 2.2m, so their detection 
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Figure 5. H band spectrum of F10214+4724. The 
position of the [Om] doublet is marked, as are 
the predicted positions of other emission lines at 
z = 2.286. 

should be less prone to difficulties in sky subtraction. 
However, it is not clear how imperfect sky subtrac- 
tion could lead to an underestimate of the H/3 flux 
without a corresponding null detection of the [Om] 
4959, 5007A doublet, since the sky lines are if any- 
thing stronger near the [Om] lines, and the relative 
variations of the sky lines are typically < 10% (Ram- 
say et al. 1992); there are also no strong atmospheric 
absorption lines in the vicinity, and the spatial varia- 
tion of the sky spectrum is also expected to be negli- 
gible. Nevertheless, the lower resolution H band spec- 
trum of Soifer et al. (1995) has a marginally detected 
blue wing to the [Om] lines which may be H/3. 

Weak features near the predicted position of the 
Hell4686 line appears to be present, which might con- 
tribute to the unusual width of the apparent 4686 
line in Soifer et al. (1995). However, none of these 
features is reliably detected, and none lies on any of 
the positions of z — 2.286 emi ssio n lines expected to 

below) which casts 



be prominent (see e.g. section 4.4 



doubt on the reliability of the apparent detections. 
Soifer et al. (1995) quote a detection of Hell 4686A at 
an Oni/Hell ratio of 23, much fainter than the limit 
(~ 5) placed in figure 5; this detection is confirmed, 
albeit tentatively, by Iwamuro et al. (1995). 

The Ha line profile in our CGS4 K band spec- 
trum (figure 6) is clearly resolved into a broader re- 
solved component, with width comparable to the UV 
emission lines, and a narrower unresolved component 
of Ha and [Nil]. We fit the K band spectrum us- 
ing the stsdas. fitting. ngauss package with the back- 
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Figure 1. WHT spectrum of F10214 + 4724. The spectrum was extracted from a full width, zero intensity aper- 
ture, thereby ensuring accurate spectrophotometry. This also included some light from a nearby companion object, 
particularly redward of 7500A, as discussed by Serjeant et al. (1995). The emission lines are marked with their 
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Figure 4. Expanded UV spectrum of F10214+4724. Positions of major emission lines are marked; also shown are 
the predicted positions of Ly/3 associated with the observed Lya. 



ground level as a free parameter, and with the fol- 
lowing four gaussians: first, three unresolved (FWHM 
200 km s _1 ) gaussians to model the narrower Ha line 
and satellite [Nil] 6548, 6584A lines, with the 6548 
component constrained to have the same redshift and 
1/3 the flux of the 6584; second, a further gaussian 
with unconstrained amplitude, position and width. 
Although the broader component is a blend of the 
broader Ha and satellite [Nil] lines, we model it as a 
single Gaussian for simplicity. 



The results, displayed in figure 7, were very simi- 
lar to a deblending with 3 Gaussians of unconstrained 
position and width, and further deblending mod- 
els with alternative simplifying assumptions were at- 
tempted. In particular, the region between the nar- 
rower Ha and the narrower [Nil] 6584 lines was very 
poorly fit in models without a broad excess, as might 
be expected from figure 7. The flux in this region is 
clear evidence for an additional component. From the 
variations in the narrower line flux between the mod- 
els we estimate the flux errors in both the narrower 
Ha and [Nll]6584A lines to be < 30% (table 1). The 



flux from the broader Ha component (distinct from 
[Nil]) is very poorly constrained. 

Note that although the narrower Ha line is unre- 
solved (FWHM less than ~ 200 km s _1 at a resolving 
power of 1428), there is marginal evidence that the 
[Nil] 6584A line is resolved. We found that at least 
some of the apparent excess on the [Nil] line, if real, 
may be attributable to the broader [Nil] component. 
If we assume that the satellite lines contribute 60% 
of this flux, we obtain a broader Ha flux of 75% of 
the total Ha. Note that the model does not imply 
the narrowest lines are Gaussian, but rather that the 
K-band spectrum must have at least one unresolved 
and at least one resolved component. This is exactly 
as expected for a hybrid starburst / AGN system. 

To summarise, we find three distinct kinematic 
components in our new spectra of F10214+4724: a 
~ 1000 km s _1 component, present in Ha and in the 
UV lines which comes from the Seyfert-II Narrow- 
Line Region as will be discussed further in Section 4; 
a < 200 km s _1 ) component, present in Ha, which we 
chose to associate with starburst activity, as will be 
discussed in Section 5; and a 4000 km s _1 Cm]1909 
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Table 1. Measurements from the WHT spectrum of F10214 + 4724. 



Line 


Position 


Rcdshift 


EW 


Line flux 


Accuracy 


Line width 


Comments 


Line 


/A 




/A 


/10- 19 W m" 2 




kms -1 




O VI 1031.9 


3388 


2.283 


60 


4.5 


30% 


1200-1600 




O VI 1037.6 


3408 


2.284 


60 


4.3 


30% 


500-1200 




CIII? 1175.7 


3866 


2.288 


6 


0.4 


60% 






Lya 1215.7 


3987 


2.279 


45 


5.3 


30% 


0-900 


two Gaussian 


Lya 1215.7 


4002 


2.292 


40 


4.7 


30% 


0-600 


fit poor 


N V 1240.1 


4074 


2.285 


130.0 


14 


20% 


1500-1700 


doublet marginally resolved 


CII 1335.3 


4390 


2.288 


4 


0.3 


60% 






Si IV 1393.7 


4577 


2.284 


8 


0.7 


50% 






O IV] 1402.5 


4607 


2.285 


40 


3.3 


30% 


1300-1500 


blended with Si IV 1402.8 


N IV] 1486.5 


4879 


2.282 


45 


4.0 


25% 


1250-1450 




C IV 1549.0 


5087 


2.284 


150 


19.0 


20% 


900-1200 




[NelV] 1602.0 


5257 


2.282 


10 


1.1 


40% 


1000-1200 




Hell 1640.5 


5385 


2.283 


100 


10.0 


20% 


800-1150 




N III] 1750 


5743 


2.282 


20 


2.0 


30% 


1200-1350 




Mg VI 1806 


5924 


2.280 


15 


1.5 


40% 


1200-1350 




Al III 1857.4? 


6104 


2.286 


5 


0.5 


60% 




Unidentified red wing 


C III] 1908.7 narrow 


6267 


2.283 


60 


6.5 


30% 


800-1000 




C III] 1908.7 broad 


6244 


2.271 


35 


5.3 


40% 


Ri3700 




[Na V] 2068 


6791 


2.265 


10 


1.2 


40% 


1000-1150 


see text 


N II] 2142.8 


7045 


2.287 


8 


0.8 


50% 


950-1100 


blend with Si VII 2148 ? 


C II] 2326.3 


7638 


2.284 


30 


3.0 


40% 


«1500 


blend with OIV 2321 ? 


[Ne IV] 2422.0 


7955 


2.284 


40 


5.2 


30% 


900-1050 




[O II] 2470.3 


8115 


2.285 


12 


0.8 


30% 


1000-1150 




[O III] 4959 


16279 


2.283 


280 


23.4 


30% 


< 1300 




[O III] 5007 


16453 


2.286 


580 


49.8 


30% 


< 1300 




[N II] 6584 


21634 


2.286 


14 


6.6 


30% 


<200 




Ha 


21564 


2.286 


12 


5.5 


30% 


<200 




Ha+[N II] broad 


21585 


~ 2.28 


83 


38.3 


20% 


~ 1800 


not dcblcnded 



Notes to Table 1: Errors on the line fluxes represent ~90 % confidence intervals expressed as a percentage of the 
best line flux estimate; for the strongest lines these are dominated by roughly equal contributions from uncertainties 
in fixing the local continuum level, and from the absolute flux calibration. Line widths were estimated from the 
FWHM of the best Gaussian fit to each line; for the rest-frame UV lines, the range lower value assumes that the 
line-emitting region fills the 2 arcsec WHT slit, the higher value that it is broadened only by the seeing. 



line which originates in a hidden Seyfert-I, or quasar, 
nucleus as has been discussed by Goodrich et al. 
(1996). 



4 THE SEYFERT-II NARROW-LINE 
REGION 

4.1 Resonant scattering of Lyman series 
lines 

There is strong evidence for resonant scattering of the 
Lya emission line: the profile of Lya is double peaked, 
and the peaks are symmetrically displaced about the 
predicted position of Lya at z = 2.286 with equal 
flux (within the errors; table 1). If the Lya-emitting 
region is free of dust then the velocity shift v of either 
Lya peak from 1215. 7Ais given by: 



v = 195( 



10 24 m- 2 



)l/3 ( 



W 4 K 



) 1/6 km B" 



(1) 



(Neufeld & McKee 1988), where N H is the total 
atomic hydrogen column density to the source of 
the Lya photons, and T is the gas temperature. 
The measured velocity separation of the peaks is 
2x ~ 560 km s _1 yielding N H = 2.5 x 10 25 m -2 for 
T = 10 4 K. 

The total line-centre optical depth to the source 
of the Lya photons is TL yQ1 ~ Nh x k-^ ya where fcLya is 
the absorption cross-section to the line centre of Lya. 
In the general case of photon absorption promoting 
an electron upwards in energy from level i to level j, 
the line-centre cross-section ki is given by 

1 V^e 2 fij 



ki = 



(2) 



47reo mcAi/D 

where fij is the oscillator strength for the transition, 
rn and e are the electronic mass and charge, c is 
the speed of light and eo is the vacuum permittivity; 
Ai/ D = (^f^) ' 5 x v is the Doppler width of the 
line with vo the central frequency, k the Boltzmann 
constant and Ma the atomic mass. In the case of the 
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2.13 2.14 2.15 2.16 2.17 2.18 2.19 
Wavelength (microns) 



Figure 6. K band spectrum of F10214+4724. The 
position of the narrower Ha and [Nil] 6584A lines 
are marked assuming z = 2.286, as is the predicted 
position of the [No] 6548A line. Note the clear 
broad base to these lines. 




2.14 2.15 2.16 2.17 2.18 

Wavelength (microns) 



Figure 7. Fit to the K band spectrum of 
F10214+4724 discussed in the text. 



Lya line, taking T = 10 4 K, these equations yield 
fc L yo = 5.8 x 1CT 18 m 2 , and thus r LyQ = 1.5 x 10 s . 

Resonant scattering is extremely efficient at re- 
moving Ly/3 photons even at moderate optical depths 
(e.g. Netzer 1975, Osterbrock 1989), in agreement 
with the null detection of Ly/3 in F10214+4724. The 
strong influence of resonant scattering on the ob- 
served properties of the Lyman series lines is a ma- 
jor difference between the spectroscopic properties 
of F10214+4724 and NGC 1068. Ultraviolet spec- 
troscopy of the latter (Kriss et al. 1992) shows a 
strong single-peaked Lya line as well as Ly/3; the ob- 
served LyP I Lya flux ratio of 0.12 implies TL yQ ~ 10. 

Is this neutral material within the narrow line 
region, or external, such as a damped Lya system 
associated with the host galaxy? The QSO proxim- 
ity effect (e.g. Giallongo et al. 1996) suggests that 
neutral material is often destroyed within the ionisa- 
tion cone. If the Nh is associated with the host, it 
need not cause damped absorption in F10214 itself 
(countering the most obvious criticism), because the 
geometry would allow scattering into the line of sight. 
This neutral hydrogen could be intrinsic to the host, 
or could be the result of accretion of gas-rich dwarfs; 
neutral systems with column densities > 10 25 m~ 2 
are known to undergo a strong cosmological evolu- 
tion which dominates the evolution in f2 g (Wolfe et 
al. 1995). 

Can the Nh gas:dust ratio, derived from the 
Lya strength, be used to resolve this? Lya photons 
are conserved in resonant scattering, but the Lya 
strength can be supressed by dust absorption between 
scatterings. The predicted case B Lya flux from our 
broader Ha flux implies a supression factor of < 10 2 , 
the limit assuming no [Nil] 6548, 84. In the damped 
Lya models of Chariot & Fall (1991), a neutral screen 
with velocity dispersion 10 km s _1 and column den- 
sity as derived above yields a Lya suppression of 
~ 10 4 , but if the gas:dust ratio is ~ 4 times lower 
than the values typically inferred from the reddening 
of quasars with damped systems this can this be rec- 
onciled with our observed Lya. The Lya flux there- 
fore cannot exclude the possibility of a damped sys- 
tem associated with the host. Alternatively, the same 
models predict a supression ~ 10 _1 for sources dis- 
tributed throughout the Nh (for a velocity dispersion 
of 1000 km s" 1 the supression becomes only ~ 0.75), 
and the authors argue that geometrical factors could 
reduce this to unity. 

There are difficulties with resonant scattering 
within the narrow line region: the clouds in the nearer 
ionisation cone suffer resonant scattering, but we 
nevertheless ought to have an unobscured view of 
the clouds in the more distant ionisation cone. The 
contradicts e.g. the observed lack of Ly/3, but this 
could be explained by differential magnification of the 
nearer ionisation cone, or selective obscuration of the 
more distant ionisation cone e.g. by the host galaxy. 
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We nevertheless favour models with at least some 
of the neutral column in the rear of the narrow line 
clouds themse lves , since the photoionisation models 
below (section 4.2) favour ionisation-bounded clouds. 
Higher resolution K-band spectroscopy may resolve 
the AGN Ha from the satellite [Nil] lines, and the 
inferred Lya supression may exclude the damped 
system model. Spatially resolved spectroscopy of 
F10214+4724 (while subject to the uncertainties in 
the projection to the source plane) may provide a fur- 
ther confirmation: diffuse Lya could be used to trace 
any extended neutral hydrogen (Villar-Martin et al. 
1996), modulo the uncertainties in the projection to 
the source plane. The detection of non-resonance lines 
spatially coincident with any extended Lya would 
rule out resonance scattering "mirrors" and identify 
such emission with the AGN narrow line region or 
starburst knots. 



4.2 Photoionisation models 

If we exclude the resonantly scattered Lya, the 
marginally resolved doublet Nvl240, 1243, and the 
blended emission lines (labelled in table 1), then we 
find no correlation of ionisation potential with emis- 
sion line width. As a result we may attempt to model 
the emission line spectrum with constant density sin- 
gle slab photoionisation models, at least to first order. 

We modelled the Seyfert-H lines with the pho- 
toionisation code CLOUDY (version 84.12, Ferland 
1993), assuming for simplicity that the rest frame 
UV emission lines in F10214+4724 suffer zero ex- 
tinction (section i.4). Unless otherwise stated, we 



used an AGN ionising continuum as characterised by 
Matthews & Ferland (1987) with a break at 10/im (see 
Ferland 1993). We considered hydrogen densities 8 < 
log 10 n < 16, where n is the total hydrogen column 
density (i.e. molecular, neutral and ionised) in units 
m -3 , and ionisation parameters —4 < log 10 U < 0, 
where U is the dimensionless ratio of incident ionising 
photons to hydrogen density, i.eU — <&(H)(nc). Here 
$ is the surface flux of ionising photons in m — 2 s — l , 
and c the speed of light. The calculation was stopped 
at a column density of 10 27 m -2 , or at a temperature 
of 4000K, since below this temperature the emission 
line flux is negligable. We used the default solar metal- 
licity (see Ferland 1993), and CLOUDY assumes a 
plane-parallel geometry. 

The Hell:Lya:Ovi:Clv ratios are sensitive to the 
UV-soft Xray continuum (e.g. Krolik & Kallman 
1988), so are sensitive to the presence or absence 
of an active nucleus. A black body ionising contin- 
uum at 40, 000K, resembling the continuum from the 
hottest OB stars, as severe difficulty in producing 
Ovil032, 1037A or Mgvi 1806A, due to the extremely 
high ionisation potentials of 4+ and Mg 4+ (.11 and 
.14 keV respectively). The AGN ionising continuum 
defined above had no difficulty in reproducing these 
lines. We are therefore confident that the ~ 1000 



km s _1 emission line region is predominantly pho- 
toionised by the active nucleus. 

None of the single-slab models provided an ad- 
equate fit to all the rest-frame UV emission lines. 
We found that similar, but not identical, conditions 
(log 10 n ~ 10, log 10 U ~ —1.5) were implied by emis- 
sion lines [Neiv] 1602, Hen 1640, Niv] 1486, Nm] 
1750, Nev 3426, Nem 3869. For example, the largest 
[Neiv] 2424 / Hen 1640 ratio in our models was with 
a density and ionisation parameter of log 10 = 10, 
log 10 U = —1.4. This underpredicts the Ovi and 
Mgvi lines, as well as the N lines. However, the Nm] 
1486, 1750 and Nvl240 ratios are all consistent with 
log 10 U ~ —1.25; this self-consistency of the N lines 
suggests that the strength of the N 1240A line may 
be due mainly to an abundance effect, as also argued 
in the IRAS detected Seyfert II NGC 3393 (Diaz et 
al. 1988) rather than due to differential magnification 
(Broadhurst & Lehar 1995). If so, this would further 
reduce the proportion of Ha in the broader compo- 
nent of figure 7. 

This model reproduced the Civ:Ciii] ratio, but 
overpredicted their strength by ~ x5. A low carbon 
abundance (~ 0.5 x solar) might resolve the anomaly, 
perhaps caused by depletion onto dust grains. Much 
more problematic is the lack of an [Om] 1665A emis- 
sion line, which in these conditions should have a flux 
~ 10% of Hen 1640A. There does not appear to be 
any simple solution; we note however that an iden- 
tical problem was found by Kriss et al. (1992b) in 
NGC1068, who argued that shock-heated gas would 
also produce the [Om] 1665 line. 

Finally, we can use our limits on the density and 
ionisation parameter to derive the sizes and masses of 
the NLR clouds. The clouds are ionisation bounded 
in the conditions we derived from the emission line 
ratios, as is also conventionally assumed for AGN nar- 
row line regions, and unlike e.g. the matter-bounded 
extended emission line region models of Wilson et 
al. 1997. This sup port s our model for the resonant 
scattering (section |4.l[ ), and is also consistent with 
the presence of [Ol] in the Keck spectrum (Soifer et 
al. 1995). Assuming the absorbing column does in- 
deed occur within the NLR clouds (assumed constant 
density for simplicity), we obtain a characteristic size 
scale of 



S = Nh/uh = 2.5x10" 



pc 



(3) 



where the errors represent the acceptable range in the 
models, rather than e.g. gaussian noise. The mean 
mass of the NLR clouds is given by 



M 



u x riH x S 



1 x 10" 



Mp 



(4) 



Unfortunately it is not possible to derive a cov- 
ering factor or volume filling factor, without spa- 
tially resolved spectroscopy. Briefly, N clouds with 
covering factor C at mean distance D parsecs from 
a luminosity L quasar, are indistinguishable from 
JV, C/4, 2D, 4L. We can however estimate the num- 
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ber of narrow line clouds, and the total cloud mass 
containted in the narrow line region. Using CLOUDY, 
we find the emissivity of a single cloud of density 



10 1 



m with ionisation parameter 10 



in e.g. the Hell 1640 line is e = i(r 2 ' 4±0 - 7 W m' 2 . 
For a magnification factor of IO.M10, our observed 
line flux implies the number of narrow line clouds is 



N = (10~ 18 W m" 2 x 47rD£)/(47reS' 2 x lQMio) 
= 1 x 10 7±1 ' 4 Alro^50 2 

where Dl is the luminosity distance to F10214+4724, 
given (assuming Q = 1, A = 0, H — 50h 50 km s _1 
Mpc" 1 ) by 



Dr 



0.04c/i" 1 (l - (l + 2.286)~ a5 )(l + 2.286) (6) 



(We also note that variations in S correlate with vari- 
ations in the cloud emissivity, so the "errors" do not 
simply add). The inferred total mass of ionised gas 
within the narrow line region is consistent with esti- 
mates in local AGN (e.g. Osterbrock 1993): 

NM = 10- ia A-KD 2 L un H S 3 /(ATvS 2 £WMio) 

= Nh x 5 x 10 12±07 m 2 M^h' 2 kg (7) 
= 2 x ltf ±Q - 7 M^KoM Q 

Our choice of Hen, while free of metallicity effects, 
assumes negligable extinction (section 4.4 ) which will 
be testable with higher quality infrared spectra. Note 
that if spatially resolved spectroscopy detects the 
emission line counterimage, it may be possible to de- 
termine TMio independently of spatial structure in the 
IRAS galaxy. As a corrolary to the covering factor cal- 
culation one may then also obtain a robust estimate 
of the luminosity of the obscured quasar (Goodrich et 
al. 1996). 

4.3 Optically-thick Ovi emission 

The resonant Ovi 1034, Nv 1240 and Civ 1549 lines 
are doublets in the isoelectronic lithium sequence. If 
such lines emanate from a region which is optically- 
thin to resonant scattering, then their doublet line 
ratios are given by the ratio of the statistical weights 
of the upper levels, i.e. 2:1 with the lower wave- 
length line of the doublet the brighter. However, in 
optically-thick regions, photon trapping and collision 
de-excitation can cause the levels to thermalize so 
that the flux ratio approaches ~1:1; this thermaliza- 
tion occurs only at high values of optical depth and 
electron density n e , namely when n e ro i5 10 22 m -3 
(Hamann et al. 1995). A doublet ratio of 1:1 there- 
fore implies extreme densities and/or optical depths. 

The Ovi 1031.9, 1037.6 doublet in F10214+4724 
is clearly resolved (Fig. 4), and its flux ratio 1:1 sug- 
gests line thermalization. We first consider whether 
this ratio is strongly influenced by blanketing by the 
Lya forest lines — both a simple empirical test and 
a brief analytical argument suggest it is not. The test 
involved fitting low-order polynomials to the UV con- 
tinuum longward of Lya, and extrapolating these fits 



(5) 



to shorter wavelengths to estimate the level of line 
blanketing. We then superimposed a model 2:1 dou- 
blet at several positions on this absorption spectrum 
and found it impossible to obtain an apparent 1:1 
ratio. This is in good agreement with the following 
analytical argument. The number of Lya forest lines 
per unit redshift per unit rest-frame equivalent width 
is well fitted by an expression of the form 



_™L - A (i + Z )i expf-— ) 

dzdw ~ w* (y + z) P( W' } 



(8) 



where 7 ~ 1.89, W ~ 0.27A and A ~ 10 (e.g. 
Bechtold 1994). By integrating this over the redshift 
range 1.770 to 1.792 (i.e. in the neighbourhood of the 
doublet and exactly encompasssing the width of either 
line of the doublet), and further integrating dN/dW, 
it can be shown that a maximum of two Lya forest 
lines are expected to lie on the position of one Ovi 
line. In order for a 2:1 doublet to be suppressed to 1:1, 
the line blanketing over the shorter wavelength line 
must be approximately a factor of 2 (in addition to 
the mean blanketing) ; the probability of this occuring 
is negligibly small, a conclusion which is insensitive to 
the uncertainties in 7, W* and Aq. 

This seems to imply that the Ovi emitting region 
in F10214+4724 is both optically thick (rovi > 1) and 
dense enough for the doublet to have become ther- 



malised, 



I.e. 7l TOvi 



10- 



We will next show 



in the that the optical depth rovi is unlikely to ex- 
ceed 10 6 , resulting in apparent broad-line-region-like 
conditions for the narrow Ovi emitting region. 

The latest release of CLOUDY (version 90.03a, 
Ferland 1996) incorporates the Ovi 1032,1037A dou- 
blet lines separately, so we used this code to explore 
the doublet ratios and Ovi optical depths for a wide 
range of physical conditions. We used the same AGN 
ionising continuum, metallicity and stopping condi- 
tions as above. The results are shown in fig. 8. We 
find that the 1:1 ratio, also where the n e ro ~ 10 22 
m -3 limit applies, is possible only with hydrogen den- 
sities in excess of riH — 10 17 m~ 3 . This conclusion 
was found to apply equally to Nn = 10 26 m~ 2 col- 
umn and/or Z — WZq metallicity gas. Such a high 
value of riH is characteristic of the broad-line region of 
AGN (e.g. Ferland et al. 1992) rather than the lower 
values (n e ~ 10 11 m~ 3 ) found in classical narrow-line 
regions. 

Such a radical conclusion can be avoided, albeit 
with some fine tuning, by considering the Ly/3 reso- 
nant absorption. This absorption is likely to be satu- 
rated, so the Ovi doublet may lie inside its damping 
wings. Using the Nh value derived above, we obtain 
an optical depth of about 0.3 for the 1032A line and 
0.09 for the 1038A line (i.e. a factor ~ 3.6). To ob- 
tain a 1 : 1 ratio from 2:1, one needs an optical depth 
of around 1 to the 1032A line, close to the estimate 
from the observed Nh- This explanation avoids the 
high densities, but is sensitive to the assumed neutral 
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Contours: n e r 0VI (in 3 ) 




Density n H (ni 3 ) 



Figure 8. Results of the photoionisation model discussed in the text. The contours show the product of the n p - 
wcightcd average electron density and Ovi 1032+1037A optical depth, n c X rovi, as functions of ionisation parameter 
and density. (Almost identical results were obtained replacing this n e with the innermost zone electron density, the 
outermost, or the total Hydrogen density.) The contours are spaced logarithmically in steps of 0.5. Only where 
n c X tqvi ~ 10 22 m~ 3 does the Ovi doublet form a 1:1 ratio. 



column, and furthermore takes no account of geomet- 
rical effects (Vilar-Martin et al. 1996). 

Inspection of the ultraviolet spectrum of 
NGC1068 (Kriss et al. 1992) suggests that the Ovi 
doublet is again in a ~1:1 flux ratio. However, there 
is no evidence of resonant scattering in this object. 



and/or Anv4720A lines can contribute up to three 
times the flux of the Hen4686A line in densities of 
10 10 — 10 12 m~ 3 , not unreasonable for AGN narrow 
line regions. Moreover, the Soifer et al. H— band spec- 
trum has very low spectral resolution, and the Hell 
4686A line flux must be treated with some caution. 



4.4 Reddening of the Seyfert-II nucleus 

We can place limits on the extinction in the observed 
UV from the Hen 1640A/Hen 1086 A ratio, which in 
the absence of reddening is predicted to be 7 (e.g. 
Seaton 1978). Although the signal to noise short- 
ward of Lya in Fig. 4 is poor, and despite possi- 
ble blanketing by Lya forest lines, the Hen 1086A 
line is marginally detected at (~7±50%)~ 1 times 
the flux of the 1640A counterpart. If real, this sug- 
gests an extinction of between zero and 0.75 magni- 
tudes around 1200A in the rest frame of F10214+4724 
(i.e. A v < 0.2), or 2600A in that of the lens (i.e. 
Ay < 0.4). We can also obtain a rough extinction es- 
timate from the Ha:Henl640A ratio, predicted to be 
about 2.3 in our solar metallicity models above. We 
find (Ha+[Nn]):Heii= 3.8, consistent with an Ha:Nn 
ratio resembling e.g. NGC 1068 (e.g. Bland-Hawthorn 
et al. 1991) and zero extinction. 

Comparison of the H band detection of the Hen 
4686A line (e.g. Soifer et al. 1995) with the 1640A 
line suggests extinctions Ay ~ 1.3. However, our 
photoionisation models above found the Fen4658A 



This low reddening in the rest-frame UV is in 
marked contrast to the high (Ay > 5) rest frame op- 
tical extinction derived from Balmer decrements (e.g. 
Soifer et al. 1995), which led several authors to infer 
that at least two physically distinct regions contribute 
to the observed emission line spectrum (e.g. Elston et 
al. 1994, Soifer et al. 1995). The limit on the Balmer 
decrement of Ha:H/3 > 20 reported by Elston et al. 
(1994) is at least in part resolved by our A"— band 
spectrum. The narrower Ha component contributes 
~ 25% of the total Ha flux,; furthermore our total 
Ha flux is ~ 25% lower than that reported by Elston 
et al., though this is within their stated photometric 
errors. However, this still yields Av > 4 in the frame 
of F10214+4724, or A v > 6 in that of the lens. (Any 
< 200 km s _1 H/3 leads to even higher A v .) Alter- 
natively, if we assume the Iwamuro et al. (1995) H/3 
detection is correct, the data are consistent with zero 
extinction. Further near-infrared spectroscopy is re- 
quired to resolve this issue. We conclude that as yet 
there is no conclusive evidence for any significant red- 
dening of the Seyfert-II nucleus. 
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5 STARBURST PROPERTIES 

Our limit on the width of the narrowest Ha compo- 
nent in F10214+4724 is suggestive of a star-forming 
region as has previously been suggested by Kroker et 
al. (1996). Moreover, the fluxes of the [Ol] 6300 A and 
[Sll] 6724A lines (Soifer et al. 1995) relative to the nar- 
rowest Ha component are similar to both those seen 
in spectra of local starburst galaxies {e.g. De Robertis 

6 Shaw 1988) and in high luminosity IRAS galaxies 
{e.g. Leech et al. 1989). We can obtain an estimate of 
the star formation rate R from the narrower Ha line 
flux S, using the following expression adapted from 
Kennicutt (1983): 

We assume a Hubble constant of Ho = 50 /150 km s _1 
Mpc -1 with Q,o = 1, Oa = 0. If we assume the star- 
burst region undergoes a net gravitational magnifica- 
tion of M. ~ 10 (Downes et al. 1995, Green & Rowan- 
Robinson 1996, Graham et al. 1995), and that aper- 
ture corrections are negligable, we obtain a star for- 
mation rate of ~ 2§h~^ 2 Mq per year. This corresponds 
to a starburst luminosity of > 2.6 x 10 11 h^ 2 Lq {e.g. 
Moorwood 1996), i.e. ~ 0.5% of the magnification- 
corrected bolometric power output. 

Reddening by intervening dust would further in- 
crease this esimate of the star formation rate. Ra- 
diative transfer models of F10214+4724 (Green & 
Rowan-Robinson 1996) suggest a much larger star- 
burst bolometric fraction of ~ 0.2 — 0.4, with a high 
starburst UV optical depth of t\jv ~ 800, i.e. an 
Ay ~ 160. Assuming this dust is well mixed with the 
Ha emitting gas {e.g. Thronson et al. 1990), the ex- 
tinction and star formation rate is quite consistent 
with our observed narrow Ha flux, which for this 
Ay yields a starburst bolometric fraction of ~ 0.5. 
In both NGC1068 (Green & Rowan-Robinbson 1996) 
and F10214+4724, the starburst and active nuclei ap- 
pear to make comparable contributions to the bolo- 
metric power output. 



6 CONCLUDING REMARKS 

The presence of the optically thick Ovi 1032, 1037A 
doublet appears at face value to imply extremely high 
densities for narrow emission line gas, tih ~ 10 17 m~ 3 . 
However, we argue that it is more easily attributable 
to the Ly/3 damping wings of resonant scattering ma- 
terial, probably withing the AGN narrow line region. 
Differential magnification appears to play a significant 
role in the emission line spectrum of F10214+4724, 
in which we identify three distinct kinematic compo- 
nents: quasar broad lines (~ 4000 km s _1 ), Seyfert II 
narrow lines (~ 1000 km s _1 ) and a starburst com- 
ponent (;$ 200 km s^ 1 ). The density, ionisation pa- 
rameter, number and total mass of Seyfert II narrow 



line clouds all resemble local Seyferts, and our inter- 
pretation of the resonant scattering agrees with that 
of Villar-Martin et al. (1996) for high-redshift radio- 
galaxies. The flux from the narrowest Ha component 
is in excellent agreement with radiative transfer mod- 
els which comprise similar quasar and starburst bolo- 
metric contributions. 
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7 APPENDIX A: COMPANION 
GALAXIES 

Figure Al shows the spectrum of a galaxy ~ 19" 
south-west of the IRAS galaxy, which fortuitously 
lay on the slit in our WHT observations. The fea- 
ture at around 5510A is a probably a low energy cos- 
mic ray event, but the emission lines at 7155A and 
(marginally) at 5323A appear to be real, since they 
occur in the spectra at both nights. Cross-correlating 
this spectrum with old galaxy models (a 1-Gyr burst 
aged by 6 x 10 s years to 1 Gyr) from Bruzual & Char- 
lot (1993) yields two redshift estimates, 0.428 and 
0.476 (both ± ~ 5%) We prefer the former, which 
identifies the weak emission features as [On] 3727 and 
[Om]5007. Close et al. (1995) report a similar red- 
shift, 0.429 ± 0.002 for a further galaxy ~ 23" south- 
west of the IRAS galaxy. The UV upturn in figure Al 
may indicate a recent starburst, as expected perhaps 
if the galaxies are tidally interacting. However, there 
could be a slight systematic error in the redshift of 
the Close et al. companion. The emission line wave- 
lengths in their spectrum of F10214+4724 are offset 
by a factor ~ 0.99 compared to other published spec- 
tra (e.g. Goodrich et al. 1995, Rowan- Robinson et al. 
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Figure A2. The sum of the spectra of sources 2 and 3 (in the notation of Matthews et al. 1994). These spectra were 
extracted with 3 pixel (~ 1.1") apertures, differing slightly to the method of Serjeant et al. 1995, though with very 
little quantitative difference. Source 2 is corrected for contamination by the IRAS galaxy by subtracting a spectrum 
of F10214+4724 scaled by the strength of the Cm] 1909 line. The sum has been smoothed with a 5-pixcl (~ 14A) 
boxcar. The flux has units 1(T 20 W m~ 2 A" 1 . 



1993), and the residuals from the 5577A sky feature 
also appear to be offset by a similar amount. 




4000 5000 6000 7000 

Wavelength (angstroms) 



Figure Al. Serendipitous galaxy ~ 19" south-west of 
F10214+4724. The flux has units 10~ 20 W m~ 2 A" 1 , 
and the spectrum is extracted with a full width zero 
intensity aperture. 



et al. 1995) is consistent with both its velocity dis- 
persion (e.g. Broadhurst & Lehar 1995) and surface 
brightness profile (Eisenhardt et al. 1996). Further 
support for this lens redshift comes from figure A2, 
in which the spectra of sources 2 and 3 (Serjeant et 
al. 1995) are summed. Source 2 is the lensing galaxy 
and source 3 is the companion ~ 3.3" north east of 
the IRAS galaxy, in the notation of Matthews et al. 
(1994). If the galaxies are indeed associated, as sug- 
gested by their similar spectral energy distributions 
and tentative 4000A breaks (Serjeant et al. 1995), 
as well as hints of tidal interaction in HST imaging 
(Eisenhardt et al. 1995), then the summed spectrum 
should show a ~ v2 improvement in signal-to-noise, 
which is indeed the case in figure A2. 

Several groups have also noted the slight overden- 
sity of galaxies in the vicinity of F10214+4724. The 
discovery of a Mgn absorber in the Keck spectrum 
at z = 1.316 (Goodrich et al. 1995), and the possibly 
associated pair at z ~ 0.4, implies that several phys- 
ically independent systems, at a variety of redshifts, 
contribute to the overdensity in the field. 



The tentative lens redshift of z ~ 0.90 (Serjeant 
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8 APPENDIX B: THE MYSTERY LINE 
AT 2067A 

In this section we summarise our attempts to iden- 
tify the line at 2067Ain the rest-frame spectrum of 
F10214+4724. We have made an extensive literature 
search, and investigated synthetic spectra generated 
by the CLOUDY photoionisation code over wide ranges 
of ionisation parameter and abundances (v84.12, Fer- 



land 1995; and see section 4.2). This led to only two 



possible identifications for this line. The first of these 
is the resonant boron Bill 2068 doublet, the next 
resonance line in the Lithium isoelectronic sequence 
below the bright OVI, NV and CIV lines. The tiny 
cosmic abundance of boron means that this line is 
predicted to be immeasurably faint, and we there- 
fore reject this possibility. The second possibility is 
the forbidden [NaV] 2067.9/2069.8 doublet (Mendoza 
1983) - a line which is not found in the CLOUDY v84.12 
output. 

To determine whether the [NaV] 2067.9/2069.8 
doublet is a plausible identification for the mystery 
line we have compared its predicted strength with 
that of [Ne IV] 2422.5/2425.1. These two pairs of 
collisionally-excited lines arise from electronic tran- 
sitions within adjacent ions in the Nitrogen isoelec- 
tronic sequence. In the low-density limit the ratio of 
the cooling rates in these doublets is given by 

— X2069 

-^2069 _ njyqv e k T n 2069 2 4 24 

L 242 4 ~~ n New X e ^24 x n 242 4 X 20 69 ( ' 

(Osterbrock 1989), where xa is the exitation po- 
tential of the upper energy level(s) and Q\ is the ef- 
fective collision strength (see Mendoza 1983). 

Adopting T = 10 4 K and solar abundances, 
and assuming that the fraction of Na in the Na 4+ 
state is similar to the fraction of Ne in the Ne 3+ 
state, the value of this ratio is m 0.005. The mea- 
sured ratio is s=s 0.23 so at first sight the puta- 
tive [Nav] 2067.9/2069.8 doublet would need to be 
anomalously strong. However, at higher densities, col- 
lisional de-excitation may become important. The 
critical densities for collisional de-excitation 

TLcrit are 

~ 10 11 m~ 3 and ~ 10 12 m -3 for the lower wave- 
length lines in the Ne and Na doublets respectively 
(in each case the higher wavelength doublet line has 
an order of magnitude lower n cr u than its part- 
ner). Given the weakness of the [OHJ3727 line, with 
ricrit ~ 10 10 m~ 3 , in F10214+4724, and the strengths 
of the [OIII] 4959/5007 line and forbidden neon lines 
(with n C rit ~ 10 12 m~ 3 ) (Soifer et al. 1995) it is cer- 
tainly feasible that collisional de-excitation can ac- 
count for the strength of the putative [Nav] doublet 
with respect to the [Neiv] doublet. 

We conclude that the [Nav]2067.9/2069.8 dou- 
blet is a probable identification for the mystery line 
in the spectrum of F10214+4724, and, by analogy, in 
NGC1068 and other active galaxies. 
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